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ABSTRACT 

We investigate the properties of both the prompt and X-ray afterglows of gamma-ray 
bursts (GRBs) in the burst frame with a sample of 33 Swift GRBs. Assuming that 
the steep decay segment in the canonical X-ray afterglow lightcurves is due to the 
curvature effect, we fit the lightcurves with a broken power-law to derive the zero time 
of the last emission epoch of the prompt emission (ti ) and the beginning as well as the 
end time of the shallow decay segment (^2 and t^). We show that both the isotropic 
peak gamma-ray luminosity (-Zjpeak,7) and gamma-ray energy (Eiso^^) are correlated 
with the isotropic X-ray energy (i?iso,x) of the shallow decay phase and the isotropic 
X-ray luminosity at t2 (I/x.ta)- We infer the properties of the progenitor stars based 
on a model proposed by Kumar et al. who suggested that both the prompt gamma- 
rays and the X-ray afterglows are due to the accretions of different layers of materials 
of the GRB progenitor star by a central black hole (BH). We find that most of the 
derived masses of the core layers are Mc = 0.1 ~ SA/q, and their average accretion 
rates in the prompt gamma-ray phase are Mc — 0.01 ~ IM0/S, with a radius of 
rc = 10® ~ 10^° cm. The rotation parameter is correlated with the burst duration, 
being consistent with the expectation of coUapsar models. The estimated radii and 
the masses of the fall-back materials for the envelope layers are re — 10^" ^ 10^^ cm 
and Me = 10~^ ^ 1-Mq, respectively. The average accretion rates in the shallow decay 
phase are correlated with those in the prompt gamma-ray phase, but they are much 
lower, i.e.. Me — 10~® ^ The values are smaller than the photospheric 

radii of Wolf-Rayet (WR) stars. In our calculation, we assume a uniform mass of the 
central BH (A/bh — IOM0). Therefore, we may compare our results with simulation 
results. It is interesting that the assembled mass density profile for the bursts in 
our sample is well consistent with the simulation for a pre-supernova star with mass 
M = 25M0. 
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1 INTRODUCTION 

One of the unexpected findings with the X-Ray Telescope 
(XRT) on-board the gamma-ray burst (GRB) mission Swift 
is the discovery of a canonical X-ray lightcurve, which shows 
successively four power-law decay segments with superim- 
posed erratic flares (Zhang et al. 2006; Nousek et al. 2006) . It 
starts with an initial steep decay following the prompt emis- 
sion. This phase usually lasts hundreds of seconds and could 
be generally explained as the tail emission of the prompt 
GRB due to the curvature effect (Kumar & Panaitescu 2000; 
Qin et al. 2004; Qin 2008; Zhang et al. 2006; Liang et al. 
2006; Zhang et al. 2007a, 2009). A shallow decay segment, 
which lasts from hundreds to thousands of seconds, is usu- 



ally seen following the GRB tail (O'Brien et al. 2006; Liang 
et al. 2007). It transits to a normal decay segment or a sharp 
drop (Troja et al. 2007; Liang et al. 2007). 

Phenomenologically, the canonical lightcurves are well 
fitted with a two-component model (Willingale et al. 2007; 
Ghisellini 2008), but the physics that shapes the canonical 
XRT lightcurves is unclear (Zhang 2007). The origin of the 
shallow decay segment is under debate. The normal decay 
segments following the shallow decay segment are roughly 
consistent with the forward shock models (Willingale et al. 
2007; Liang et al. 2007), favoring the long lasting energy in- 
jection models for the shallow-decay segments (Zhang 2007). 
The chromatic transition time observed in both the X-ray 
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and optical afterglows challenges this scenario (Fan & Piran 
2006; Panaitescu 2007; Liang et al. 2007). 

Alternative models were proposed to explain the shal- 
low decay segment (see review by Zhang 2007). The shal- 
low decay would result in a high gamma-ray efficiency (e.g. 
Zhang et al. 2007), and loka et al. (2006) proposed that 
the efficiency crisis may be avoided if a weak relativistic 
explosion occurs 10'' — 10^ s prior to the main burst or 
if the microphysical parameter of the electron energy in- 
creases during the shallow decay. Shao & Dai (2007) in- 
terpreted the X-ray lightcurve as due to dust scattering of 
some prompt X-rays(c.f., Shen et al. 2009). The scattering 
of the external forward-shock or of the internal shock syn- 
chrotron emission by a relativistic outflow could also explain 
the observed X-ray afterglows (Shen et al. 2006; Panaitscu 
2007). Uhm & Beloborodov (2007) and Genet, Daigne & 
Mochkovitch (2007) interpreted both X-ray and optical af- 
terglow as emission from a long-lived reverse shock. Liang 
et al. (2007) argued that the physical origin of the shallow 
decay segment may be diverse and those shallow decay seg- 
ments following an abrupt cutoff might be of internal origin 
(see also Troja et al. 2007). Ghisellini et al. (2007) suggested 
that the shallow-to-normal transition in the X-ray afterglows 
may be produced by late internal shocks, and the transition 
is due to the jet effect in the prompt ejecta (see also Nava 
et al. 2007). Racusin et al. (2009) also suggested that the 
shallow-to-normal transition may be a jet break occurring 
during energy injection. Interestingly, Yamazaki (2009) re- 
cently suggested that the X-ray emission might be an inde- 
pendent component prior to the GRB trigger. By shifting 
the zero time point of the shallow-to-normal decay segment 
in the canonical XRT lightcurves, Liang et al. (2009) found 
that the shallow-to-normal decay behavior might be due to 
a reference time effect. 

It has long been speculated that long GRBs are associ- 
ated with the deaths of massive stars and hence supernovae 
(SNe) (Colgate 1974; Woosley 1993; see Zhang & Meszaros 
2004; Piran 2005; Meszaros 2006; Woosley & Bloom 2006 
for reviews) . The coUapsar model is the most promising sce- 
nario, in which the GRB jets are powered by the accretion 
of an accretion disk or a torus fed by the fall-back mate- 
rial from the coUapsar envelope (e.g., Popham et al. 1999; 
Narayan et al. 2001; Kohri & Mineshige 2002; Di Matteo et 
al. 2002; Kohri et al. 2005; Lee et al. 2005; Gu et al. 2006; 
Chen & Beloborodov 2007; Liu et al. 2007; Kawanaka & Mi- 
neshige 2007; Janiuk et al. 2007; Janiuk & Proga 2008). Ku- 
mar et al. (2008a, b) proposed that the canonical lightcurves 
may be produced by the mass-accretion of different layers 
of progenitor stars. In the framework of their model, the X- 
ray emission of GRBs may give insight into the properties 
of the progenitors. In this paper we investigate the charac- 
teristics of the X-ray afterglow lightcurves in the GRB rest 
frame and infer the properties of progenitor stars with a 
sample of 33 GRBs based on the model of Kumar et al. Our 
sample selection and the method are presented in §2. In §3, 
we give the correlations between the prompt gamma-rays 
and the X-rays in the shallow decay segment. Inferred pa- 
rameters of progenitor stars are reported in §4. The results 
are summarized in §5 with some discussion. Throughout, a 
concordance cosmology with parameters Hq = 71 km s~^ 
Mpc"\ = 0.30, and = 0.70 are adopted. 



2 DATA 

The XRT data are downloaded from the Swift data archive. 
The HEAsoft packages, including Xspec, Xselect, Ximage, 
and Swift data analysis tools, are used for the data re- 
duction. We use an IDL code developed by Zhang et al. 
(2007b) to automatically process the XRT data for all the 
bursts detected by Swift/BAT with redshift measurements 
up to October of 2008. Our sample includes only those XRT 
lightcurves that have a clear initial steep decay segment, a 
shallow decay segment, and a normal decay segment. We 
get a sample of 33 GRBs. We fit the spectra accumulated 
in the steep and shallow segments with an absorbed power 
law model and derived their spectral indicefl Regarding the 
steep decay segment as a GRB tail due to the curvature ef- 
fect (e.g., Liang et al. 2006, Zhang et al. 2007a, 2009; Qin 
2009), we estimate the time of last emission episode of the 
GRB phase with the relation a = 2-\- (5 (Kumar et al. 2000; 
Liang et al. 2006). We fit the steep-to-shallow decay segment 
with, 

F = Fo[(^)-(2+''^' + (!l^)-(2+/3i) X {-y^-l (1) 

tl tl t2 

where ti is zero time point of the last emission epoch of 
the prompt gamma-rays and t2 is the starting time of the 
shallow decay segment. The end time of the shallow decay 
segment (ts) is taken as the break time between the shallow 
to normal decay phases. Flares in the steep-to-shallow decay 
segments are removed, if any. Technically, the shallow decay 
segments are poorly sampled for some GRBs. We flx the 
Q2 value in order to get a reasonable flt. Illustrations of 
our fitting results for 24 bursts of our sample are shown 
in Figure. 1. We derive the X-ray fiuence Sx in the time 
interval [t2, ia] in the XRT band and calculate the isotropic 
X-ray energy with £iso,x ~ A-kD\Sx /{^ + z), where Dl 
is the luminosity distance. The isotropic peak fluxes of the 
prompt gamma-rays (Lpoak,7) are in 1024 ms timescale. We 
take the X-ray luminosity at t2, Lx,t2, as a characteristic 
luminosity of the shallow decay segment. Our results are 
summarized in Table 1. With the data reported in Table 1, 
we show the distributions of ti, t2, and t^, in comparison 
with GRB duration Tgo in Figure 2. It is found that the 
distribution of t\ is comparable to Tgo, t2 is about 100—1000 
seconds, and t^ is in lO"* - 10^ seconds. 



3 CORRELATIONS 

The correlations between the prompt gamma-rays and the 
X-rays in the shallow decay segment may reveal some phys- 
ical relations between these two phases. We show the pair 
correlation of the observables between the two phases in 
Figure 3, and measure these correlations with the Spear- 
man correlation analysis. Our results are reported in Table 
2. We flnd that there are several outliers at Tgo < 30s in 
the correlation of Tgo and ti. It seems natural that for long 
bursts tl could be a mark of the end of the prompt emission 
epoch and will be very likely approximately equal to Tgo, 

^ Although the steep decay segment has significant spectral evo- 
lution (Zhang et al. 2007b), we derive only the time-integrated 
spectral index for our analysis. 
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since the time of the last pulse should occur close to the end 
of the overall emission. This breaks down for shorter bursts, 
where the offset between Tgo and ti becomes important. We 
also find that Tgo is not correlated with the time intervals 
t2 — ti and ts — t2, indicating that the break features in 
the XRT lightcurves are independent of the durations of the 
prompt gamma-rays. The energy releases in the two phases 
are strongly correlated. The Z/peak,7 and -£150,7 ^re correlated 
with the X-ray luminosity at t2 (Lx,t2)- This fact likely sug- 
gests that the two phases may have related energy budgets 
from the same central engine, and the physical conditions to 
power the gamma-rays and the X-rays should be similar. 



4 PROPERTIES OF PROGENITOR STARS 

Kumar et al. (2008a, b) proposed that the prompt gamma- 
rays and X-ray afterglows are due to accretion of difi'erent 
layers of a coUapsar by a newly-formed black hole (BIf) with 
IOM0 . In their model, the highly variable lightcurves of the 
prompt gamma-rays are explained as production of the ac- 
cretion of the dense, clumpy materials of the stellar core, 
and the power-law decay X-rays may be due to the accre- 
tion of the fall-back materials of the progenitor envelope. In 
this section, we derive the properties of the progenitor stars 
based on the model of Kumar et al. 

As mentioned by Kumar et al. (2008b) and from the 
numerical simulations from GRMHD by McKinney (2005), 
the efficiency of the accreted energy to the radiation is likely 
to depend on many details. Here we assume a uniform radi- 
ation efficiency of 1% of the accreted mass by a rotating BH 
(McKinney 2005). We also do not consider the beaming ef- 
fect. Then the masses accreted by the BH during the prompt 
gamma-ray phase (M-y) and during the shallow decay phase 
(Mx) are estimated with Mace ~ 100i?iso/c^. The average 
accretion rate thus can be estimated with M ~ Macc/Tacc, 
where Tacc is the accretion timescale in the rest frame. Con- 
sidering the fall-back of total accreted particles as free-fall, 
the radius r for the fall-back time T can be estimated with 



rio 



(2) 



where no = r/lO^" cm, Mbh,i = Mbh/WMq, and = 
T /lO^s. We assume Mbh,i = 1 in this work. The rotation 
rate fn (r) of the fall-back material at radius r is defined as a 
ratio of the local angular velocity fl{r) to the local Keplerian 
velocity f2fc(r) of the material at r. 



Mr) 



nfc(r) 



(3) 



Considering the viscosity among accreted particles before 
they reach the BH and combining with Eq. (2), one can 
obtain 



fn{r) (X ( 



lOT' 



\l/3 



(4) 



where Ovia is the viscous parameter and face is the viscous ac- 
cretion timescale of the fall-back material. Please note that 
the timescale Tacc is different from face as tacc ~ 2/avisf2fc is 
the viscous accretion time of the fall-back material after it 
has circularized but Tacc is the accretion time for fall-back 
material within difi'erent layers of the progenitor star, which 



is given by the fall-back time without considering the viscos- 
ity among the particles. Assuming that the observed fiux is 
proportional to the accretion rate and that the timescale t^u 
of the decrease of mass fall-back rate from /2 to /i is much 
larger than tacc, the upper limit of f(i{r) can be obtained by 
(Kumar et al. (2008a, b) 



f-acc Jl 



(5) 



A lower limit on fnir) is derived from the centrifu- 
gally supporting condition that the fall-back material is able 



to form an accretion disk at a radius ra 
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Td > S-Rg, where Rg = GMbh/c^ (c is the velocity of the 
light). For a convective envelope, the density profile at r is 
p oc r~*, where 5 is determined by the slope of the shallow 
decay segment, with 5 = 8(0:2 -I- l)/2. The sharp decline in 
the steep decay indicates that the density in the transition 
region decreases sharply. In this region, tacc 3> tdu, so the 
accretion in this region can be ignored. 

With the data of the bursts in our sample, we calculate 
the radii for stellar core Tc, transient region rt, and envelope 
region the limits of spin parameter /n, the index of the 
density profile (5, the accreted masses {Mc and Me), and the 
average accretion rates (Mc and Me) in the prompt gamma- 
ray and the shallow decay phases. They are tabulated in 
Table 3. We show the distributions of these parameters in 
Figure 4. We find that the derived radii of the core layers of 
the progenitor stars for all the bursts are Tc = 10^ ~ lO'^" 
cm with the rotation parameter as /n,c = 0.02 ~ 0.05. The 
masses of the core layers for about two-thirds of GRBs in 
our sample are Mc = 0.1 ~ 5M0 with a mass density of 



10^ 



10 g cm , and their average accretion rates in the 



prompt gamma-ray phase are Mc — 0.01 ~ IMq/s. 

For the envelope layer, the estimated radii, lower lim- 
its on the rotation parameters, and the masses of the fall- 
back materials are = 10^° ~ 10^^ cm, the lower limit 



for the rotation parameter /o^e = 10 



10" 



and Me = 



10^ ~ IM0, respectively. The average accretion rates in 
the shallow decay phase are much lower than that in the 
prompt gamma-ray phase, i.e.. Me = 10~* ~ lO~^M0/s, 
but they are correlated. We measure the correlation with 
the Spearman correlation analysis, which yields log Mc = 
(-4.03 ±0.21) -I- (0.78 ±0.14) log Me with a correlation coef- 
ficient r = 0.72 and a chance probability p < 10~*, as shown 
in Figure 5. The estimated mass density in the envelope is 
g cm^ 



5 CONCLUSIONS AND DISCUSSION 

We have investigated the characteristics of the X-rays in the 
GRB rest frame, and inferred the properties of progenitor 
stars with a sample of 33 GRBs based on the model of Ku- 
mar et al. Assuming that the steep decay segment is due 
to the curvature effect, we fit the lightcurves with a broken 
power-law to derive the zero time of the last emission epoch 
of the prompt emission (ti), the beginning {12) and the end 
time (ts) of the shallow decay segment. The Tgo is roughly 
consistent with ti, but it is not correlated with the time 
intervals of t2 — ti and tz ~ t2. The -Eiso,7 and Lpcak,7 are 
correlated with Eiso,x and Lx,t2 ■ This fact likely suggests 
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that the energy budgets for the two phases may be from the 
same central engine. 

Based on a model proposed by Kumar et al. (2008a, b), 
we inferred the properties of the progenitor star with both 
the prompt gamma-rays and the X-ray data. The derived 
radii of the core layers of the progenitor stars for all the 
bursts are Vc = 10^ ~ 10^° cm with a rotation parameter as 
fn,c = 0.02 ~ 0.05. The masses of the core layers for about 
two-thirds of GRBs in our sample are Mc = 0.1 ~ 5Mq 
with a mass density of 10^ ~ 10^ g cm , and their av- 
erage accretion rates in the prompt gamma-ray phase are 
Mc = 0.01 ~ IMq/s. The estimated radii, lower limits on 
the rotation parameters, and the masses of the fall-back ma- 
terials for the envelope layers arc re = 10^" ~ 10^^ cm, 

= 10"^ ~ 10"^ and Me = 10"^ ~ IMq, respectively. 
The average accretion rates in the shallow decay phase are 
much lower than those in the prompt gamma-ray phase, i.e., 
Me = 10"* ~ W^Mq/s, but they arc correlated. The esti- 
mated mass density in the envelope is ~ 10^* g cm~^. 

The connection between long-duration GRBs and SNe 
was predicted theoretically (Colgate 1974; Woosley 1993) 
and has been verified observationally through detecting 
spectroscopic features of the underlying SNc in some nearby 
GRBs (Woosley & Bloom 2006). The collapsar model is the 
most promising scenario to explain the huge release of energy 
associated with long duration GRBs (Woosley & Weaver 
1995; Paczyhski 1998; MacFadyen & Woosley 1999; Zhang 
et al. 2003; .Tauiuk & Proga 2008). In this scenario GRBs arc 
produced by a jet powered by accretion of the core and the 
fall-back materials of the progenitor star through a torus. 
We infer the properties of the progenitor stars by assuming 
that both the prompt gamma-rays and the X-rays observed 
with XRT are due to the accretion of different layers of pro- 
genitor stars. We compare the distribution of re, the radius 
of envelope region, with the photospheric radii of a sample 
with 25 WC-type and 61 WN-type Wolf-Rayct (WR) stars 
(Koesterke & Hamann 1995, Li 2007) in Figure 6. We find 
that the photospheric radius of a WR star is larger than re, 
consistent with our prediction. 

In our calculation, we regard that all GRBs are from 
a unified collapsar with a central BH of Al = lOAf0. Our 
results thus might be compared with simulation results for a 
collapsar with a given mass. We compare the derived mass 
density profile as a function of radius r with simulations 
in Figure 7, in which the simulated mass density profile is 
taken from Woosley & Weaver (1995) for a massive star with 
M = 25Mq (see also Janiuk & Proga 2008). It is found 
that, although the derived p are systematically larger than 
the results of simulations, they are very consistent with the 
sirrmlations. 

In the collapsar models, the accretion duration should 
be as long as the material fall-back timescale from the col- 
lapsar envelope available to fuel the accretion disk or torus. 
Rotation of the progenitor star should be high enough to 
form the disk or torus. One thus might expect a relation 
between the burst duration and the rotation parameter fn- 
For a progenitor star with higher rotation, the angular mo- 
mentum loss should be longer, and the accretion timescale 
might be longer, hence a longer GRB event. We show the 
correlation between ti and fn in Figure 8. A tentative cor- 
relation is found, with a correlation coefficient of 0.72 and a 
chance probability p < 10""*. This correlation indicates that 



the higher fn, the longer GRB could be observed, consistent 
with the model's expectation. 

It is believed that GRBs are highly coUimated, with a 
beaming factor ft ~ 1/500 from the optical afterglow obser- 
vations(Frail et al. 2001). The measurement of the beaming 
angle has also proven exceedingly difficult in the Swift era 
(Cenko et al. 2009). In our analysis, we do not considering 
the beaming effect. The lack of detection of jet-like breaks 
in the late XRT lightcurve might suggest that the X-ray 
jet would be less collimatcd than expected from the optical 
data (Burrows & Racusin 2006; Liang et al. 2008). In spite 
of this, the accretion rates and the accreted masses could be 
up to 2 orders of magnitude lower than those derived in this 
analysis if beaming is considered. 

We appreciate valuable suggestions and comments from 
the anonymous referee. We also thank Bing Zhang, Li-Xin 
Li, Zhuo Li, and Tong Liu for helpful discussion. 
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Table 1. XRT Observations and our fitting results 
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The X-ray fluence integrated from t2 to is and its error in the XRT band (0.3-10 keV) 
^ The time-averaged photon index of the steep decay phase 

The zero time of the emission epoch corresponding to the steep decay segment 

t2 and ta are the begin and the end times of the shallow decay segment 
^ The slopes of the shallow decay segment 

Table 2. Spearman correlation coefficients between the prompt gamma-rays and the X-rays in the shallow decay phase 



Tao 0.55(8. 1E-4;|^ 0.58(2.9E-4) 0.08(0.63) 0.27(0.13) 0.46(6. 3E-3) 
-t'peak,7 0.40(0.02) 0.08(0.66) -0.30(0.08) 0.74(<lE-4) 0.66(<lE-4) 

Eiso,7 0.51(2.2E-3) 0.33(0.06) -0.18(0.31) 0.68(<lE-4) 0.72(<lE-4) 



" In the bracket is the chance probability. 
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Table 3. Inferred properties of the progenitor stars for the bursts in our sample 
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Radii of the core, transit and envelope region of the progenitor stars. 
^ The range of the rotational parameter of the core layer, /n = ^P- . 

The lower limit on the rotational parameter of the envelope layer. 

The density profile in the stellar envelope: p(r) ~ r"*. 
^ The fall-back mass during the prompt phase. 
f The fall-back mass during the shallow decay phase. 
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Fig.l (continued) 

Figure 1. Illustration of the XRT light curves (dots with error bars) with our best fit red solid for some bursts in our sample. The 
and degrees of freedom of the fits are also marked in each panel. 
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Fig.l (continued) 
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Figure 2. Distributions of logarithmic T90, ti, t2, and 43. 
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Figure 3. Correlation between the properties of the prompt gamma-ray phase and the shallow decay phase. Lines are the best fitting 
results. 
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Figure 5. Correlation of the average accretion rates between the prompt gamma-rays Mc and the shallow decay phases Me- The solid 
line is the best fit to the data with a Spearman correlation coefficient r = 0.72 and chance probability p < 10~*. 
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Figure 6. Comparison of tiie radii of the envelope layers of the bursts in our sample {step dash line ) with the photosphcrical radii of 
86 WR stars [step solid line) 




Figure 7. Assembled mass density as a function of radius for the bursts in our sample with comparison of the simulation for a pre- 
supernova star with mass of 25Mq (the curve). 
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Figure 8. Correlation between the rotational parameter and the burst duration. 



